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Abstract We follow the premise that most intermediate luminosity optical transients 
(ILOTs) are powered by rapid mass accretion onto a main sequence star, and study the 
effects of jets launched by an accretion disk. The disk is formed due to large specific 
angular momentum of the accreted mass. The two opposite jets might expel some of the 
mass from the reservoir of gas that feeds the disk, and therefore reduces and shortens 
the mass accretion process. We argue that by this process ILOTs limit their luminosity 
and might even shut themselves off in this negative jet feedback mechanism (JLM). The 
group of ILOTs is a new member of a large family of astrophysical objects whose activity 
is regulated by the operation of the JLM. 

Key words: stars: jets — (stars:) binaries: general — stars: mass-loss — accretion, ac¬ 
cretion disks 


1 INTRODUCTION 


supernovae and above those of novae (e.g. Mould et al. 

1990; Rau et al. 2007^ Ofek et al. 2008; 

Prieto et al. 2009; Botticella et al. 2009; 

Smith et al. 2009; 

Berger et al. 2009bf Kulkarni & Kasliwal 

2009 

; Mason et al. 201(1 Pastorello et al 

2010; Kasliwal et al. 2011; Tvlenda et al. 2013; Kasliwal 

2013 

) has been growing in recent years 

Kasliwal 2013). Excluding low luminosity supernovae and 


similar objects, such as Ca-rich transients and .la SNe, the rest of the eruptive events are part of a group 
of the so called intermediate luminosity optical transients (ILOTs: lBerger et aPl2009bl ). 

ILOT events of systems that harbor asymptotic giant branch (AGB) or extreme-AGB (ExAGB) pre- 
outbur st stars, like N GC 300 OT20Q8-1 (NG C 300QT: iMonard120081: iBond et alJl2009t iBereer et al.l 
2009b ) and SN 2008S dArbour & Bo lesl2008l) . were studied both in t he frame of single s t ar models (e.g.. 


Thompson et al. 2009; Kochanekll201 ll) and binary stellar models dKashi et alJ l2010fc iKashi & Soketl 


vived these events (Adams et al. 12015l) . 


2010b : Soker & Kashj 201 1[ 201 2 . 2013 ). Recent observations cast doubt that the projenitor star 


Me lev & So kef d2014i) conclude that single-star models for 


ILOTs of evolved giant stars encounter severe difficulties, and that ILOTs are most likely powered 
by a binary interaction. We here continue with developing and exploring the binary model for all ILOTs 
as we outlined in an ea r lier pap er dKashi & Sokeill2010bh . 

In IKashi & Soketl ((2010b) we developed the High-Accretion-Powered ILOT (hereafter, HAPI) 
model. The HAPI model suggests that many (most) ILOTs are powered by a high-accretion rate event 
onto a main sequence (MS), or slightly involved off the MS, star in a binary system. In some cases 
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the binary system is actually a triple star system, where the tertiary star induces orbital instabilities 
and causes the two inner stars to interact. In this mass transfer event the accreted mass possesses 
high specific a ngular momentum to form an accretion disk, or an accretion belt, around the MS star 
(IKashi & Sokeni2009l) . The accretion disk can launch two opposite jets that expel more mass from the 
system to form an expanding bipolar nebula dKashi & Soken l2010a). suc h as the bipolar nebula of 77 
Car, the Homunculus, that was formed in the Great Eruption (GE; e.g., iHumphrevs & Martini 12012 
and references therein) . The connection between some ILOTs and b ipolar planetary nebulae (PNe) was 
mentioned in th e past (I Soker & Kashi|[2oT~2t lAkashi & Sokerll 2 Q 1 3l) in regards to the binary model for 
ExAGB-ILOTs. lPrieto et all (2009) already made a connection between NGC 300OT and pre-PNe. The 
connection between some bipolar PNe and and some ILOTs has ga in ed support in recent years (e.g., 
iBoumis & Meaburnll2013t IClvne et aDl20l4 iDe Marcoll20 1 4l2015allbt IZiilstrall2015l) . 

In this paper we update the classification and grouping of ILOTs (section 0 to further emphasize 
the common properties they share, including the presence of jets. In section [3] we turn to propose that in 
some ILOTs the jets lead to a negative jet feedback mechanism (JEM). We compare some properties of 
the JEM in ILOTs to those in other systems where the JFM operates. Our short summary is in section[4] 


2 ILOTS 


2.1 The Energy-Time Diagram (ETD) 

The objects studied here (gap objects that are not supernovae) can be classified in the following way. 
ILRT: Intermediate-Luminous Red Transients. These are events resulting from evolved stars, such as 
AGB and ExAGB stars, and similar objects, e.g., red giant branch (RGB) stars. Most likely a companion 
accretes mass and the gravitational energy powers the eruption, e.g., NGC 300 OT. A similar process 
occurs in giant eruptions of luminous blue variables (LBVs). 

LBV giant eruptions and SN Impostors. Giant eruptions of Luminous Blue Variables. Examples in¬ 
clude the 1837-1856 GE of 77 Car and the pre-explosion eruptions of SN 2009ip. ILRTs are the low 
mass relatives of LBV giant eruptions. 

LRN or RT: Luminous Red Novae or Red Transients or Merger-bursts. These are powered by a 
full merger of two stars. The process of destruction of the less dense star onto the denser star releases 
gravitational energy that powers the transient. Examples include V838 Mon and V1309 Sco. Merger 
events of stars with sub-stellar objects are also included. 

ILOT: Intermediate-Luminosity Optical Transients. ILOT is the term for the combined three groups 
listed above (ILRT, LRN and LBV giant eruptions). These events, we argue, share many common physi¬ 
cal processes, in particular being powered by gravitational energy released in a high-accretion rate event, 
the HAPI model. 

We emphasize again that the definition of ILOT does not include low luminosity supernovae and 
similar objects, such as Ca-rich transients and .la SNe (for these see iKasliwalli2013h . 

A panoramic way to examine and compare ILOTs and other transient events is by using the Energy 
Time Diagram (ETD), as presented in Fig. Q] (For a diagram showing peak luminosity rather than total 
energy against the duration of the event, see, e.g., iKulkarni et al1i2007 an d lKasliwal2013l ) The diagram 
shows the total energy of the transient as a function of its duration. Both the energy and the duration are 
not easy to define. The total energy includes both the integrated luminosity (radiated energy), and the 
kinetic energy that went to the ejected material. The radiated energy is calculated by integration of the 
the light curve in the optical bands. The bolometric radiated energy is then calculated using a bolometric 
correction that is not always accurate. The kinetic energy is not easy to derive from observations, and 
may be a few times larger than the radiated energy. It results in a range of possible energies, that is 
presented by an error bar in the ETD. 

In some cases the energy that is missed by observations, i.e., the total available energy minus that 
emitted in the optical band, can be large. This is the case for V838 Mon, where the energy required to 
inflate the envelope, marked by the blue filled circle in Fig. [T] is more than 10 times larger than that 
carried by radiation dTvlenda & Sokerl2006l) . To account for the total available energy we calculate and 
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Energy-Time Diagram (ETD) 



1.5 2 2.5 

log(timescale/day) 


Fig-1 Observed transient events on the Energy-Time Diagram (ETD). Blue empty circles 
represent the total (radiated plus kinetic) energy of the observed transients as a function of the 
duration t of their eruptions. The Optical Transient Stripe (OTS), is an approximately con¬ 
stant luminosity region in the ETD. It is populated by ILOTs of different kinds: Intermediate- 
Luminous Red Transients (ILRT), Luminous Red Novae (LRN) and LBV giant eruptions, aka 
SN impostors. It also includes predicted BD-planets merger-bursts dBear et al.l201 ll) . The up¬ 
per limit of the OTS is constrained by equation ©. The green line represents nova models 
computed using lu minosity and duration from ldella Valle & Livio dl995h. Nova mod els from 
lYaron et all d2005l) are marked with red crosses, and models from Shara et ahl (1201 0t> are rep¬ 
resented with diamonds. The total energy does not include the energy which is deposited in 
lifting the envelope that does not escape from the star. Lor ILOTs that have model to estimate 
the gravitational energy released by the accreted mass (the available energy), it is marked by 
a black asterisk. PNe that might have been created in ILOT events are marked with gray hor¬ 
izontal lines, indicating the uncertainty in their ILOT event durations dSoker & Kashiff20T2l) . 


present the gravitational energy released by the accreted mass, as described in eq. (5) of lKashi & Sokeil 
d2010bi) . This energy is marked by a black asterisk in the ETD. The advantage of doing so is that the total 
available energy better reflects the physical processes behind the transient, as according to our model 
these transient events are powered by gravitational energy. However, the total energy is estimated using 
a detailed model, and those models were derived for only a fraction of all transient events. 

The ILOT time scale is defined typically as the time from its peak luminosity until it reaches its 
pre-ILOT luminosity, before the rise. (Using a criterion like a decrease by 3 magnitudes will not change 
much the graph.) Ideally it is defined for the bolometric luminosity, and as observations are limited, the 
waveband that contains most of the energy is used. 
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The ILOTs are found to lie on a wide stripe in the ETD - the Optical Transients Stripe (OTS). The 
stripe is not just an accidental position, and has its roots in accretion physics, as we demonstrate in 
section [L2l below. Being a transient object positioned on the OTS is one of the shared characteristics of 
ILOTs, and in most cases it suggests that the transient is powered by gravitation energy, and belongs to 
the ILOT group. 

In previous papers we reported differ ent ILOTs that lay on the OTS dKashi & Sokerl 120 10b: 
iBear et aill 201 lklSoker & Kashill20lHl2012h . The ETD presented in Fig. |T]is updated with a number of 
recent ILOTs, as listed below. 


1. SN 2009ip. An LBV that suffered several eruptions before (e.g.. [ Berger et~aT]l2009 aL iMaza et al.l 
120091: IPraJke etal. 2012; Levesaueetalj2014; Mauerh anetal JI2013 : Pastorelloetal. 2012 ) its last 


one i n 2012 ( Mar gutti et al.ll2014 ; Martin etld] 2015 ), that might have been a real supernova explo¬ 

sion dSmith et all2014l) . There is also a p ossibility t hat the last 2012 er uption was an exceptionally 
strong giant eruption as well (So ker & Kashil 12013k iKashi et al.ll2013l) . The calculated location of 
the sub-eruptions in 2011 and the two in 2012 ( 2012 a and 2012 b), as well as the total three years 
activity (SN 2009ip Tot), are marked on Fig. H] It has been assumed that the total energy is twice 
the radiated energy for all eruptions of SN 2009ip. 

SN 2010mc. This super nova had a pre-explosion outburst about a month before explosion 
dOfek et ali 120131) . ISokeil d201 3bh proposed that the pre-explosion outburst was powered by ac¬ 
cretion onto a main sequence companion. 

HDE 269006). A giant eruption of an LBV. Energy is calculated based on iMehner et al.l 
Late time data is missing so the location may be somewhat more to the upper right (longer 
duration and more energy), but certainly within the OTS. 

CK Vul (Nova V ul 1670). An histori c eru ption showing three p eaks between 1670 and 1672 
dShara et alJI 19851) . lHaiduk et al.l d2013l) and iKamihski et all d2015l) suggested it might have been 
an ILOT. We used historic data to position it on the ETD. Observations of a bipolar nebula 
dHaiduk et al .1120071) sug gest the activity of je ts. 

OGLE-2002-BLG-360. rfvlenda et al.l d2013l) claim it to be an ILOT (a n LRN), 

6 . SN2007sv. This seems to be a SN impostor according to the analysis of Tartaglia e t alJ d2015l) . We 
find it to be at the center of the OTS. The total radiated energy is more than an order of magnitude 
below that of typical SNe. Since more than 7 years have passed since the outburst (ILOT event) and 
no SN explosion took place, this ILOT is not excited by oxygen burning in the core of a massive 
star. Its location on the OTS shows it js one of the short duration and low energy SN impostors. 
M31LRN 2015. A recent LRN dWilliams et alll2015l) . 

SN Hunt248. This SN Impostor experienced an eruption with 3 peaks and returned to its pre erup¬ 
tion luminosity in about 230 days dKankare et al. Il201~5l and references therein). Its position in the 
ETD is close to that of the first 2012 peak of SN 20009ip. 


2 . 


3. 


4. 


5. 


2.2 The high-accretion-powered ILOT (HAPI) model 


The most significant property ILOTs share according to the HAPI model is that the energy source is 
gravitational. A high accretion rate M a onto a MS or a slightly evolved off-MS star accounts for the 
high luminosity. We here explore some of the implications of the HAPI model. 

Let M a and R a be the mass and radius of the star accreting the mass, respectively. Star ‘ 6 ’ is the 
one that supplies the mass to the accretion; it is possibly a MS star that is completely destroyed, as in 
V838 Mon and V1309 Sco, or alternatively an evolved star in an unstable phase of evolution that loses 
a huge amount of mass, as in the GE of 77 Car. The average total gravitational power is the average 
accretion rate multiplied by the potential well of the accreting star 


Lg 


GMgMg 

Ra 


( 1 ) 


The accreted mass will plausibly form an accretion disk or an accretion belt. The accretion time ought to 
be longer than the viscous time scale for the accreted mass to lose its angular momentum. The viscous 
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timescale is 




-1 


C'/v* 

0.1 


-1 


Ra 


5 R<* 


3/2 


M a 

8 Mq 


- 1/2 


days, 


( 2 ) 


where v = a C S H is the viscosity of the disk, H is the thickness of the disk, C s is the sound speed, a 
is the disk viscosity parameter, and v$ is the Keplerian velocity. We scale M a and R a in equation (0 
to confirm with to the parameters of V838 Mon (iTVlenda et all2005l) . For these parameters, the ratio of 
viscous to Keplerian timescale is \ = ivisc/iif — 160. 

The accreted mass is determined by the details of the binary interaction process, and varies for 
different objects. We scale it by M acc = ?/„ M a . Based on the modeled systems (V838 Mon, V 1309 Sco, 
77 Car) this mass fraction is r] a < 0.1 with a large variation. The value of rj a < 0.1 can be understood 
as follows. If the MS s t ar collides with a star a nd tidally disrupts it, as in the model for V838 Mon 
dTvlenda & SokeHl2006t ISoker & Tvlendalf2006h . the destructed star is likely be less massive than the 
accretor M acc < Mj, < 0.3 M a . In another possible case an evolved star loses a huge amount of mass, 
but the accretor gains only a small fraction of the ejected mass, as in the GE of 77 Car. 

The viscous time scale gives an upper limit on the accretion rate 


m < VaMa ~ 4 (Ih-\ f ( H / R A ( c °/ v <t> 

“ t visc VO.l/ vo.l ) V o.l ) V o.l 


Rg 
5 Ra 


The maximum gravitational power is therefore 


-3/2 


Mg 
8 Mp, 


3/2 


Mq yr x . 


<T max = ^. 7.7X10- 


Ra 


(-) ( — 

Vo.l ) V 160 . 


-1 


Rg 
5 Rn 


-5/2 


Mg 

8 Mq 


5/2 


erg s 


-1 


(3) 


(4) 


where we replaced the parameters of the viscous time scale with the ratio of viscous to Keplerian time 
X- Eq. (Q} determines the upper bound on the OTS in the ETD. The upper bound might be crossed if the 
accretion efficiency 77 is higher and/or the stellar parameters of the accreting star are different. For most 
of the ILOTs the accretion efficiency is lower, hence they are located below this line, giving rise to the 
relatively large width of the OTS. The uncertainty in rj a is large and may be even above unity, but only 
in extreme cases. Therefore, we seldom expect to find objects slightly above the upper limit. 


2.3 Implications of Accretion Powered Events 

Adopting the premise that ILOTs are powered by high mass accretion rate, the HAPI model, a great deal 
can be learned about them. We bring three examples. 

The stellar masses of the two stars composing the 77 Car system were constrained using that model. 
The GE o f 77 Car was initiated by close interact ion with the eccentric companion close to periastron 
passages dKashi & Sokeil2010al) .lD aminelil d 19961) tried to support the existence of a periodical variation 
using the historical data, but claimed it cannot be done due to poor confidence on data. Our calculation 
worked in an opposite way - we used the available historical light curve to develop a model and calculate 
the stellar and orbital parameters. 

According to the model of lKashi & Sokerl d2010al) . tidal forces of the companion disrupt the LBV 
that was in an unstable state, having an Eddington factor close to 1. According to the HAPI model the 
luminosity peaks of the GE resulted from accretion onto the companion close to periastron passages 
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(iKashi & Sokeri[2010a ). Those peaks are ~ 5.1 years apart, while the present orbital period is ~ 5.54 
years, suggesting that during the GE in the 1840’s, the orbital period was shorter. Changes in orbital 
period came most likely from mass loss by the erupting system (e nlarges the peri od), and accretion 
from the more massive LBV to the companion (shortens the period). [Kashi & Sokerl d201 Oal ) found that 
the only way for the HAPI mechanism to work within the orbital changes constraints, requires the two 
stars to have a significantly larger masses than initially thought. Instead of a binary system composed 
of an LBV with a mass of M\ ss 120 Mq and a companion of M2 ~ 30 Mq, the mass es should be 
in the range M\ « 170 - 200 Mq and M 2 ~ 80 Mq. Based on the results of lFiger et al.l dl998l) a 
zero-age main sequence (ZAMS) star with an initial mass of Mzams ~ 230 Mq is required to explain 
the present luminosity of the primary of 77 Car. The massive primary star deduced here is expected to 
results from a somewhat more massive ZAMS star, compatible with Mzams ~ 230 Mq. This cannot 
be said on a present primary mass of 120 Mq. 

The accretion of the mass ejected from the primary onto the secondary formed a disk and launched 
jets. These jets had two effects: (1) They shaped the Homunculus as a bipolar nebula; (2) They may have 
cleared some of the mass in the vicinity of the companion and by that reduced the amount of available 
gas for further accretion, through the JFM discussed in section[3] 

A similar idea was applied to the LBV star P Cyg. There is no observed companion to P Cyg. 
Adopting the HAPI model and analyzing the peaks in its historical light-curve. [Kashi! (i2010t) concluded 
that a 3-6 Mq binary companion in an eccentric orbit most likely does exist. The periastron passages 
of the companion caused the LBV to lose mass, part of which was accreted by the companion, releasing 
energy. The whole process modified the orbital parameters, resulting in a present period of ~ 7 years. 
The nebula around P Cyg has bi-axial features that may have been a result of relatively weak jets 
launched from the proposed companion during the outbursts. If such a companion is not found, it is 
possible that is was swallowed by the LBV star. 

[Soker & Kashil (i2Q12h proposed that a number of bipolar PNe may have undergone ILOT events. 
Some PNe, such as NGC 6302 and the pre-PNe OH231.8+4.2, M l-92, and IRAS 22036+5306 have 
bipol ar features that were formed in a rapid mass ejection event (e.g., Mea burn et alJ2008tlSzvszka et al.l 
1201 lh . The formation of a bipolar nebula in an event lasting weeks to years is a prediction of the HAPI 
model. The interaction with a MS companion causes the AGB (or extreme AGB) to lose a huge amount 
of mass. Part of this mass is accreted by the MS companion through an accretion disk. The disk launches 
two jets that form the lobes. The process is accompanied by high luminosity that makes the event an 
ILOT. 

3 THE NEGATIVE JET FEEDBACK MECHANISM 

Consider a flow structure where a large gas reservoir supplies mass to an accretion disk that launches 
jets. If these jets are not too narrow, they can interact with the reservoir and expel and heat the mass 
residing there. This heating and/or ejection of reservoir gas reduces the mass supply rate to the accre¬ 
tion disk, hence lowering the jets’ power, and even shutting-off jet formation completely. Jets can affect 
the amb ient gas very efficiently by getting shocked to high temperature and inflating hot bubbles (e.g., 
[Hillel & Sokenl2016l and references therein). The bubble has a larger cross section to push and expel 
ambient gas, also in directions away from the direction of the jet. The inflation of the bubble is accom¬ 
panied by the formation of many vortices that mix part of the hot gas with the cooler ambient gas (e.g., 
[Hillel & Sokerll20~i~6i) . This results in an efficient heating of the ambient gas. This entire process can 
convert kinetic energy to thermal energy and then to radiation. 

This negative jet feedback mechanism (JFM) was studied for a variety of astrophysical objects in 
the literature. These systems are listed in Table 1. The table does not include planetary nebulae (PNe), 
as for them the interaction of the jets is on very large scales and it does not set a feedback. Only if the 
shaping takes place close to the bi nary sy stem, like if the nebula is expelled during a GEE, that might 
also involved an ILOT (see lSokenl2016l) . feedback occurs. It is listed here under GEE. We compare 
the size and the gravitational potential of the accreting body in each type of objects with those of the 
reservoir of gas, R a , $ a , R res , and $ res , respectively. 
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Table 1 Comparing Feedback Properties in Astrophysical Objects. 


Object type 

Accreting 

object 

Ra 

(cm) 

( km s -1 ) 2 

-^res 

(cm) 

^res 

(km s -1 ) 2 

^ a / ^ res 

Galaxy 

formation 

SMBH 

10 ii_i “ 

c 2 

-j^q^U-^3 

(300) 2 

« 10° 

Cluster 
cooling flows 

SMBH 

10 13 - 16 

c 2 

10 21 - 24 

(1000) 2 

« 10 s 

Core collapse 
supernovae 

NS 

10 6 

(100 000) 2 

10 9 

(10 000) 2 

« 100 

GEE [1] 

MS star 

10 11 

(500) 2 

10 13 

(30) 2 

« 100 

CEE [1] ’ [21 

MS star 

10 11 

(500) 2 

10 13 

(30) 2 

« 100 

ILOTs [2] 

MS star 

10 11 - 12 

(1000) 2 

10 11 - 13 

(30 - 1 000) 2 

Si 3 - 100 


Notes: 

R a and Tv stand for the typical radius of the accreting object, and the gravitational potential on its surface. R les 
and T>res stand for the typical radius of the reservoir of gas for accretion and the energy required to expel it from 
the system, c represents the speed of light. 

Abbreviations: BH: black hole; CCSN: core-collapse supernova; CEE: common envelope evolution; GEE: 
grazing-envelope evolution; NS: neutron star; SMBH: super-massive BH. 

1 We refer here only to a giant primary star and a main sequence (MS) companion. See text for the case of a NS or 
a BH companion. 

2 While in the other objects jets are a crucial ingredient in the evolution, in the CEE and and in ILOTs in some 
cases jets do not occur, or play a small role. Furthermore, not all cases the JFM operates even if jets do exist. 


ISoker et all (12013b made comparisons of the JFM operating during galaxy formation, in cooling 
flows in clusters of galaxies, during the common envelope evolution (CEE), in PNe, and in core collapse 
supernovae (CCSNe). Boosted by the results of the Chandra and XMM-Newton X-ray observatories, 
there have been many papers in the last fifteen yea r s on the operation of a JFM in galaxy formation and 
cooling flow clusters (McNa mara & Nulsenll2012t ISoker et al.ll2013l and references therein). In recent 
years there have been suggestions for the operation of a JFM in the explosion of massive stars in core 
collapse supernovae (CCSNe), and during th e CEE, The jet-feedback mechanism for the explosion of all 
CCSNe is termed the jitteri ng-jets model (IPapish & Sokerll201 ll 1201 3. l2014allhl : iGilkis & Sokeil|20l4 
120 1 Sallbl IPapish et alJl2Q16h . 


Numerical simulations encountered difficulties in expelling the envelope gas during the CEE, as 
a large fraction of the envelope stays bound to the system even after a substantial spiraling-in of the 
two s tars (e.g., ISandquist et alJ l~i~998l iLombardi et al.[|2006t iDe Marco et al.l I 20 TH iPassv et all 1201 ll 
12012b iRicker & Taamll2012l bu t see iNandez et al.ll2015 ). These difficulties and the problems with the 
acE~prescription (e.g., Sokerl 120 1 3ah lead to the suggestion that in many cases envelope ejection is 
facilitated by jets launched by the more co mpact companion dSokerll2004t iKashi & Sokerll201 llISokerl 
1201 3al 120141) . I Armitage & Li viol (l2000i) and lChevalierl (120 1 2l) studied CE ejection by jets launched from 
a neutron star (NS) companion but not as a general CE ejection process. 


An efficient JFM might even prevent the formation of a CEE. Instead, the compact companion 
spirals-in while grazing the envelope of the giant star, in what is termed the grazing envelope evolution 
(GEE; [SokeifcO 15l) . If the companion is on an eccentric orbit and grazes the envelope, or even forms a 
temporary common envelope, it can accrete at a high rate, releasing gravitational energy th at is observed 
as an ILOT. As such, the JFM can operate in ILOTs, as we propose below. We note that llvanova et al.l 
( 120131 ) mention that the ejection of the envelope through a CCE can lead to an ILOT event, but they 
mention no jets. 
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Launching jets by a MS star in the GEE and the CEE requires that MS stars can accrete at high rates 
reaching 0.001-0.1M@ yr _1 , and launch jets. It seems that MS stars can indeed accrete mass at high 
rates if the jets that are launched by the accretion disk or belt remove most of the accreted energy and 
angular momentum (IShiber et alJl2016l) . As well, even if the specific angular momentum of the accreted 
gas is below the value needed to form a Keplerian a_ccretion disk, but not by much, the accretion belt 
formed around the MS star might launch jets rtSchreier & Sokeil2016l) . An accretion belt is defined here 
as a sub-Keplerian accretion flow on the surface of the accreting body, with mass concentration closer 
to the equatorial plane. ISchreier & Sokerl d2016l) argue that jets might be blown from the polar regions 
of the accretion belt. 

The effect of the JFM can be quantified by the ratio between the gravitational potential energy of the 
mass accreted onto the surface of the compact body, and that of the ambient gas in the reservoir. If this 
ratio is large, as in galaxies and clusters of galaxies, then a small amount of accreted mass can launch 
jets that have a large impact on the ambient gas. We find that for many ILOTs this ratio is smaller than 
in the other types of objects listed in Table 1. This implies that the JFM in ILOTs is less efficient that in 
most other cases listed in Table 1. A large fraction of the mass available in the reservoir is accreted in a 
typical ILOT event. 

We note that in the case of ILOTs, jets can interact with gas residing close to the accreting star, hence 
leading to the operation of the JFM, or with gas residing further out that is not part of the reservoir. The 
interaction with gas residing further out, but in a still optically thick region, converts kinetic energy to 
thermal energy, and then radiation. Namely, the interaction can increase the ILOTs luminosity. 

We turn now to discuss several settings of ILOTs powered by accreting MS stars, and the different 
manifestation of the JFM in each one of them. The three different scenarios differ mainly in the time 
period during which the accretion process takes place, and in the reservoir of the accreted mass onto the 
MS star. 

(1) An ILOT event at a periastron passage. In this scenario the accretion phase lasts for a fraction 
of the orbit. The accreted gas originates in one side of the bloated atmosphere of a giant star. For 
these, the JFM does not destroy the entire mass reservoir, and if the MS companion survives, then the 
outburst might repeat itself. The best example is the GE of r/ Car, where there were at least two outbursts 
separated by ~ 5.2 years of orbital period. Both mass removal and mass accretion act to change both 
the eccentricity and the orbital separation, but in opposite directions. Whether the eccentricity and the 
orbital separation increase or decrease depends on the detail of the interaction, and on the mas s loss rate 
far from periastron, when accretion does not take place or has a small rate (see Kashi & Sokerl[201 Obi) . 
In the GE of p Car, according to the HAPI model, the companion accreted mass mainly at periastron 
passages. The distance between the two stars at periastron passages was a p ss '■><)()R . : . The companion 
could have entered the extended and tidally distorted envelope of the primary LBV star, and then get 
out. This may be termed GEE. Most of the accreted mass came from the equatorial plane. Some mass 
though, could have been accreted from all direction including from above and below the companion, 
at a typical distance from the secondary star of « 50-100 Rq. The radius of the secondary star is 
R a ~ 20 Rq. When the mass of the primary star is included, we find $ a /i> res « 3. This relatively 
low value and the small fraction of reservoir gas residing along the polar directions make the JFM less 
efficient than in most objects listed in Table 1. Still, we hold that the JFM can operate in giant eruptions 
of very massive stars. 

(2) The JFM in a fall back gas. In this scenario the JFM takes place a long time, relative to the 
dynamical time, after the outburst was triggered. The accreted gas is a fall back gas from the violent 
event, e.g., a merger process. This JFM is expected to operate once after one violent event. As not 
much mass falls back, the JFM here will not change much the binary parameters. The estimated ejected 
mass in the 1988 outburst of M3J_ RV was « 0.001-0.1 M 0 , and with an expansion velocity of « 
100-500 km s -1 dMould et al.lll99Cil) . Based on some similarities with V838 Mon, Soker & Tvlendal 
(2003) suggested that M31 RV eruption was caused by merger-burst event. In Kashi & Soker (2010b) 
we raised the possibility that a merging processes might be replaced with a rapid mass transfer episode 
in a binary system. We considered a small ejected mass that quickly became optically thin. Part of the 
ejected material might have not reached the escape velocity, and fell back toward the star, or the binary 
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system if the companion survived the event. The inflated envelope and the gas that did not reach the 
escape velocity fell back toward the star. Instead of a long-lived inflated envelope as in the case of V838 
Mon, we argued, an accretion disk was formed after about a year. Because of the high specific angular 
momentum in the binary system the fall back gas formed an accretion disk around the accreting star. 
The material closer to the center formed an accretion disk before the outer parts of the envelope have 
collapsed. 

Here we add that such a scenario can lead to the onset of the JFM. The disk from the fall back gas 
might have launched jets. The jets could have expelled some of the outer parts of the envelope that did 
not form a disk yet, hence lowering the amount of gas accreted to the center. This is a negative JFM 
process. It might act in general cases of merger-bursts where a large envelope is inflated, and in addition 
the gas that continues to be accreted onto the central star has sufficient specific angular momentum to 
form an accretion disk that launches two jets. Namely, even if most of the mass in a red transient (a 
merger-burst) is lost via the outer Lagrange point, e.g. JPeicha et al.i (120151) . jets might still be formed. 
This can be after the merger process as outlined above, or before the full merger when the denser star 
accretes gas from the other star via an accretion disk. 

The ratio <J> a /$ res depends on the mass distribution in the envelope that is formed after merger. 
Even if it is an extended envelope, most of the reservoir mass might reside close to the accreting star, 
resulting in $ a /$ res « 10. As well, due to the high angular momentum of a merger event, mass is 
concentrated toward the equatorial plane, hence less affected by the polar jets. Over all, the JFM is not 
extremely efficient in this setting. Still, it plays some role. 

We note that in the merger simulation of Nandez et al.l (120141) no accretion disk was formed. We do 
not claim that jets are formed in all merger-bursts, but in many cases they might form. 

(3) Spiral-out grazing envelope evolution. In this scenario the accretion phase lasts for the entire 
orbit, and for many orbits. The accreted gas originates from the outer parts of the giant’s envelope and 
near its equator. Mass removal can be significant, hence acting to increase orbital separation. But as the 
JFM requires the MS companion to accrete mass from the envelope, even if orbital separation increases, 
it will not be by much, as the companion cannot get too far from the envelope. Basically, the companion 
might end at an orbital separation of about the size of the giant primary star, « 1 AU. The removal of 
gas by jets can have an indirect influence through dynamical effects. This is known to occur for example 
during galaxy formation. A rapid removal of baryonic mass by AGN activity reduces the depth of the 
gravitational well, hence causing also the dark matter to expand (e.i_ T .. [Martizzi et ap2013l) . This further 
reduces the gravitational well, and hence acts to further lower the density of gas in the center of the 
galaxy. This in turn acts to reduce accretion rate to the central SMBH. 

Mass loss by either star in a binary system and mass transfer between them change the orbital 
parameters. The operation of the JFM in removing gas from a giant envelope (the primary) can have a 
dynamical effect. The jets remove mass above and below the trajectory of the accreting compact star 
through the giant envelope, as depicted in Fig. [2] As the center of mass is between the two stars, in the 
case that the secondary star is not much lighter than the primary star, the removed mass is closer to an 
axis through the center of mass and perpendicular to the equatorial plane (depicted by a dashed line) 
than most of the rest of the giant’s surface. This implies that the removed envelope mass has a lower 
value of specific angular momentum than the average of the giant’s surface. Since in a strong binary 
interaction it is expected that the mass loss will be concentrated to the equatorial plane, we find that 
the JFM removes mass with a lower specific angular momentum than that from an equatorial mass loss 
from the primary star. Therefore, the mass removal by jets has a less pronounced effect on reducing 
the orbital separation. In addition, the energy used to expel the envelope mass comes from the mass 
accretion onto the companion, and not from the orbital energy. These two properties, of a relatively low 
specific angular momentum of the removed gas and using the accretion energy to remove gas, do not act 
strongly to reduce the orbital separation. Therefore, the mass loss process in the JFM acts to increase 
the orbital separation. On the other hand, if the accreting secondary star is much lighter than the primary 
giant star, mass transfer acts to reduce the orbital separation. 

If the accretion rate started at a high rate and later decreased, there may be enough gas in the disk to 
lunch a jet. The mass removal by jet can supersede the decreasing accretion rate, deplete the reservoir 
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Fig. 2 A schematic flow structure of a grazing-envelope evolution (GEE). The companion 
accretes mass through an accretion disk and launches two jets in the outskirts of the giant 
envelope. The jets inflate bubbles and efficiently remove most of envelope mass residing 
above and below the secondary orbit. A large fraction of this mass is located not far from the 
axis through the center of mass, depicted by a dashed line perpendicular to the orbital plane, 
and hence posses a low value of specific angular momentum. 


and consequently stop the mass transfer. This would cause an increase of the orbital period. If the binary 
has an eccentric orbit, eccentricity can increase as well. This spiral-out evolution resulting from the JFM 
operating in binary stars will be studied in detail in a forthcoming paper. 

4 SUMMARY 

We explored some aspects of the physics of ILOTs under the assumption that they are powered by high 
mass accretion rate on to a MS, or slightly evolved off the MS star. In this high-accretion-power ILOTs 
(HAPI) model one star in a binary system accretes mass from a companion. The companion can survive 
or be completely destroyed in a merger-burst event. 

In section l27T1 we defined ILOTs as the general group of gap objects, between novae and supernovae, 
that are not type .la SNe or related objects. We also classified ILOTs into three subgroups according to 
the mass transfer process. We updated the ETD where ILOTs occupy a stripe, termed the OTS (Fig. 
[Q). The upper boundary of the OTS can be accounted for by accretion via an accretion disk around a 
MS star. This upper boundary was derived in section I2i2l (eg. l4l). Some implications of the HAPI model 
were discussed in section 12731 In particular, we reiterate an earlier result that the masses of the two stars 
composing the LBV binary system 77 Car are much more massive than what is usually assumed based 
only on the Eddington luminosity limit. 
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Our most significant new claim is that in some cases ILOTs are regulated to some degree by a 
negative feedback mechanism mediated by jets. In this JFM the jets expel some mass from the ambient 
gas reservoir of the accretion disk. This in turn reduces and/or shortens the mass accretion rate. In Table 
1 we compared some properties of the JFM of ILOTs with other types of astrophysical objects. In section 
[3]we discussed three examples of the possible operation of the JFM in ILOTs. 

The examples listed in section [3] show that the JFM in ILOTs is less efficient than in most other 
types objects listed in Table 1 for two reasons. ( i ) The value of the ratio of gravitational potential on the 
mass accreting body to that in the reservoir $ a /<b res , is typically lower than in other types of objects. 
This implies that more jets’ mass is required to remove a given amount of reservoir gas. ( ii ) Because 
the interaction occurs relatively close to the accreting objects, the high value of the specific angular 
momentum of the accreted gas implies that it tends to flow closer to the equatorial plane. As the jets are 
lunched along the polar directions, they interact with a smaller fraction of the ambient reservoir gas than 
in a spherical distribution. 

As more ILOTs are being discovered we continue to populate the ETEQand learn more about them. 
Further work is required to quantify the outcome of the JFM. Numerical simulations will be the next 
step. These can show that the jets actually casue the depletion of the mass reservoir. Most relevant is the 
study of the newly proposed GEE, that requires 3D hydrodynamical numerical simulations. 
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